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Abstract 

We present general constraints on dark matter stability in hadronic decay channels 
derived from measurements of cosmic-ray antiprotons. We analyze various hadronic 
decay modes in a model-independent manner by examining the lowest-order decays 
allowed by gauge and Lorentz invariance for scalar and fermionic dark matter 
particles and present the corresponding lower bounds on the partial decay lifetimes 
in those channels. We also investigate the complementarity between hadronic and 
gamma-ray constraints derived from searches for monochromatic lines in the sky, 
which can be produced at the quantum level if the dark matter decays into quark- 
antiquark pairs at leading order. 



1 Introduction 



Despite the overwhelming evidence for the existence of dark matter in our Universe, very 
little is known about its properties from the point of view of particle physics. One of the 
most striking properties of the particles comprising the dark matter is their longevity, 
evidenced by their survival from production in the early Universe to the present day. 
Among the hundreds of matter particles which have been produced and studied over the 
last century at accelerators, only very few have lifetimes on a cosmological scale: the 
proton, the electron and the three neutrinos. The longevity of the electron, which is the 
lightest charged particle, can be attributed to the conservation of electric charge. The 
longevity of the lightest neutrino, which is also the lightest fermion, is due to Lorentz 
symmetry. The longevity of the other two neutrinos, which can decay into the lightest 
neutrino and a photon, is because of their tiny masses. Lastly, the longevity of the proton 
can be attributed to the accidental conservation of baryon number in the renormalizable 
Standard Model Lagrangian, which might be broken by higher-dimensional operators. 
However, no general principle is known which would necessitate a cosmologically long 
dark matter lifetime. On the contrary, the longevity of the dark matter particle can only 
be understood within concrete models of dark matter, and is usually due to additional 
symmetries imposed from the start. 

In the absence of a fundamental reason to guarantee the absolute stability of the 
dark matter particle, it is important to analyze the potential signatures of the decay of 
dark matter particles and the constraints on the dark matter lifetime which follow from 
cosmic-ray measurements. Generally, indirect constraints on dark matter stability can 
be inferred from observations of gamma rays [Tj [2J El IU El E], neutrinos El El EE], 
electrons/positrons [TTl El US] and antiprotons/antideuterons [T6l H7j US} fl9l 120] . 

In this paper we focus on the hadronic decay modes of general dark matter candidates. We 
use the measurements of the antiproton-to-proton fraction performed with the PAMELA 
satellite instrument [211 E2] to constrain the minimum dark matter lifetime in hadronic 
decay modes, examining a number of possibilities for the coupling of the dark matter 
particle to the Standard Model particles. We also argue that the partial lifetime for 
the decay of dark matter particles into a quark-antiquark pair and a neutral fermion, 
V'dm QQN, can also be constrained using gamma-ray observations; more concretely via 
searches for monoenergetic gamma-rays which can be generated at the quantum level by 
radiating photons from virtual quarks. Instead of photons, weak gauge bosons can also 
be generated at next-to-leading order in this manner, even in models where the dark 
matter does not decay hadronically at leading order. 

This paper is organized as follows. In Section [2] we review the production and prop- 
agation of antiprotons in the Galaxy. We then derive constraints on the partial dark 
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matter lifetime in various hadronic decay modes in Section |3j In Section H] we discuss 
the generation of monochromatic gamma rays from decays into quarks and compare the 
resulting constraints with those from cosmic-ray antiprotons. The production of antipro- 
tons from radiative decays into weak gauge bosons is analyzed in Section We present 
our conclusions in Section |6j 

2 Antiproton Production from Dark Matter Decay 

Antiprotons are a rare component of the cosmic radiation (at 10 GeV kinetic energy 
about one antiproton is observed for every 10 4 protons [2Tj). These cosmic antiprotons are 
believed to be overwhelmingly or exclusively of secondary origin by spallation of primary 
cosmic rays on interstellar gas, a hypothesis which is in good quantitative agreement with 
observations [23]. Therefore, the cosmic-ray antiproton flux constitutes a sensitive probe 
for particle production in the Galaxy from exotic sources like dark matter decay. In this 
section we briefly review the production and propagation in the Galaxy of antiprotons 
from dark matter decay. 

Given a population of dark matter particles with mass ttidm and lifetime tdm, the 
production rate of antiprotons per unit energy and unit volume at a position f with 
respect to the center of the Galaxy is given by 



where dNp/dE is the energy spectrum of antiprotons produced in the decay and Pdm(^0 
is the density profile of the distribution of dark matter in the Milky Way halo. For 
definiteness we will adopt the spherically symmetric Navarro-Frenk- White halo density 
profile [Ml [25]: 



with po — 0.26GeV/cm 3 and r c ~ 20 kpc, which corresponds to a local dark matter 
density of p & = 0.4 GeV/cm 3 [261 123 [28] . Due to the linear dependence on the dark 
matter density, different choices for the halo profile yield very similar constraints [29] . 

After being produced in the Milky Way halo, antiprotons propagate through the 
Galaxy and its diffusive magnetic halo in a rather complicated manner before reaching 
the Earth. Cosmic-ray propagation in the Galaxy is commonly described using a sta- 
tionary two-zone diffusion model with cylindrical boundary conditions [30]. Due to the 
relatively large antiproton mass, we can neglect energy losses as well as reacceleration 
effects in our treatment [31] . We furthermore neglect catastrophic losses from antiproton 
annihilation with ordinary matter in the Galactic disk since they have a negligible effect 



Qp(E,r) 



Pdm(^) dN p 
^dmTdm dE 
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on the total antiproton flux. Under these approximations, the number density fp(T, f, t) 
of antiprotons per unit kinetic energy T satisfies the following transport equation: 

= ^ = V • [K(T, r)Vf p - V c (r)f p ] + Q P (T, r) . (2.3) 

The first term on the right-hand side of the transport equation is the diffusion term, 
which accounts for the propagation through the tangled Galactic magnetic field. The 
diffusion coefficient K(T, r) is assumed to be constant throughout the diffusion zone 
and is parametrized by K(T) = K (3TZ 5 , where (3 = v/cis the antiproton speed as a 
fraction of the speed of light and 1Z is the rigidity of the particle, which is defined as the 
momentum in GeV per unit charge, 1Z = p(GeV)/\Z\. The second term is the convection 
term, which accounts for the drift of antiprotons away from the disk induced by the 
Milky Way's Galactic wind. Following Ref. [32] we will assume that it has axial direction 
and that it is constant inside the diffusion region above and below the Galactic disk: 
V c (r) = V c sign(z) k. Lastly, Qp(T,r) is the source term of cosmic antiprotons as defined 
in Eq. ([2~11 . 

The solution of the transport equation at the heliospheric boundary at our position 
in the Galaxy, r = r , z — 0, can be formally expressed by the convolution 

Tin ax 

o 

where T max = m-DM — Tn p , with m p being the proton mass. Analytical and numerical ex- 
pressions for the Green's function Gp(T, T') derived within the semi-analytical framework 
described above can be found in 1291. 



Given the number density, the flux of primary antiprotons from dark matter decay 
at the heliospheric boundary is finally given by 

$ dm (t) = v M m , (2.5) 

where v is the speed of the antiprotons as a function of their kinetic energy. 

The computation of the local antiproton flux suffers from various sources of uncer- 
tainty. The most important one stems from the choice of the propagation parameters 
that enter in the transport equation (the normalization of the diffusion coefficient, the 
power-law index of the rigidity, the height of the diffusion zone and the speed of the 
convective wind). This uncertainty in the model parameters amounts to as much as two 
orders of magnitude. In view of this large uncertainty we calculate our constraints for the 
full range of parameters compatible with observations. In contrast, the uncertainties on 
the choice of the dark matter halo profile and the local dark matter density are subdom- 
inant, and amount at most to a factor of order one [29]. Due to the linear dependence of 
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K n (kncVMvr) 


L (kpc) 


V, (km /si 


MIN 


0.85 


0.0016 


1 


13.5 


MED 


0.70 


0.0112 


4 


12 


MAX 


0.46 


0.0765 


15 


5 



Table 1: Astrophysical parameters compatible with the observed B/C ratio that yield the 
minimum (MIN), medium (MED) and maximum (MAX) flux of antiprotons [32]. 

the fluxes on the dark matter density, our constraints can easily be rescaled to other val- 
ues. Therefore, we will just present results for the Navarro-Frenk- White halo profile with 
Po = 0.4 GeV/cm 3 , whereas we will present the entire range of results for three different 
propagation models that are consistent with the observed boron-to-carbon (B/C) ratio 
and that give the maximum (MAX), medium (MED) and minimum (MIN) antiproton 
flux [32J. The relevant parameters are summarized in Tab. [TJ 

Lastly, antiproton propagation inside the Solar System is influenced by solar modula- 
tion effects, which have an effect on the locally observed low-energy antiproton spectrum. 
Under the force field approximation [33], the antiproton flux at the top of the Earth's 
atmosphere is related to the interstellar antiproton flux [34] by the simple relation: 

®* (Tt ° a) " I 2^ + 7* ) ^ (Tls) ' (2 ' 6) 

where Tis = Ttoa + <f>F, with Tjs and Ttoa being the antiproton kinetic energies at the 
heliospheric boundary and at the top of the Earth's atmosphere, respectively, and (f)p 
being the solar modulation parameter, which varies between 500 MV and 1.3 GV over 
the eleven-year solar cycle. In order to compare our predictions with the PAMELA data 
we will take (p F = 550 MV [35]. 

3 Tree-Level Decays 

We begin our study of the constraints on the dark matter lifetime by examining the 
tree-level decays of scalar and fermionic dark matter particles into weak gauge bosons, 
Higgs bosons and quarks, which subsequently hadronize, producing antiprotons and other 
particles. 

3.1 Decay of Spin-1/2 Dark Matter Particles 

We first consider the decays of a spin-1/2 dark matter particle. In this case, the simplest 
decays allowed by Lorentz invariance are two-body decays into a spin-1/2 fermion and 
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a boson (scalar or vector) and three-body decays into fermions. Decays with larger mul- 
tiplicities in the final state, such as three-body decays into a spin-1/2 fermion and two 
scalars, four-body decays into three fermions and one gauge boson or five-body decays 
into fermions will not be considered in this paper. To keep the analysis as general as 
possible we will not specify the nature of the neutral fermion in the final state, which we 
denote by N. This particle could be a neutrino or a heavier particle, possibly another 
dark matter component, such as a gravitino or a gaugino of some hidden-sector gauge 
group. 

The lowest order decay modes of a spin-1/2 dark matter particle are: 

• Decay into a spin-1/2 fermion and a vector boson: ^dm — > Z°N, ipr>M 

• Decay into a spin-1/2 fermion and a Higgs boson: ipr>M — > h°N 

• Decay into two quarks and a spin-1/2 fermion: xjj — > UiUjN, ^dm — > didjN, ipDM — > 



Here, i,j = 1,2,3 are quark generation indices (note that the lepton generation is irrel- 
evant for antiproton production). In the following we only regard the case that the two 
quarks are of the same flavor, % — j. 

The production of antiprotons by fragmentation of the weak gauge bosons, the Higgs 
boson and the quarks was simulated using the event generator PYTHIA 6.4 [36J. Taking 
these spectra as input in the source term in Eq. 12.11 we then calculate the antiproton-to- 
proton fraction at Earth by solving the transport equation, Eq. f)2.3p . for the MIN, MED 
and MAX set of propagation parameters that we list in Table [H From the resulting 
antiproton-to-proton ratios we calculate lower bounds on the dark matter lifetime for 
dark matter masses in the range 100 GeV - 20 TeV by comparing the results to the 
PAMELA measurements of the antiproton-to-proton fraction. We have also adopted the 
background flux of secondary antiprotons produced in spallations of cosmic rays in the 
interstellar medium calculated in [23] which was calculated within the same two-zone 
diffusion model and agrees well with the measured antiproton-to-proton ratio. We derive 
conservative constraints on the partial dark matter lifetime from the requirement that 
the total antiproton-to-proton ratio should not exceed the PAMELA measurements at 
95% CL. as determined by the x 2 statistic 



where the Xi are the measured values of the antiproton-to-proton ratio at kinetic ener- 
gies Ti, and the <7j are the corresponding statistical errors. The $(Tj) denote the top-of- 
atmosphere fluxes of protons and antiprotons. The total flux $ tot (T) = $ DM (T) + $ bkg (T) 
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(3.1) 
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is the sum of contributions from dark matter decay and the "background" from ordi- 
nary astrophysical production by cosmic-ray spallation. We also regard the case of a 
pure antiproton signal from dark matter decay, where we disregard the astrophysical 
background, and for which we use the following "x 2 " statistic, 

x 2 = J2 WdM^)/^^)-^} 2 j (3. 2) 
i < ~ r% 

where R(x) = x for x > and R(x) = for x < 0, so that only fluxes that exceed 
the measurements are penalized when we regard the case of a pure dark matter-induced 
signal with no background. 

Given that astrophysical models of secondary antiproton production by cosmic-ray 
spallation can quantitatively reproduce the observed antiproton-to-proton fraction well, 
the constraints derived in the former analysis are significantly more stringent than in the 
latter. 



3.1.1 Two-Body Decays 

We analyze the decay of a spin-1/2 dark matter particle in a simple toy model where we 
assume a Lagrangian that describes the interaction between dark matter, weak gauge 
bosons and fermions, where we assume the dark matter to be a Majorana particle. The 
two-body decay of a spin-1/2 dark matter particle with mass wdm into a weak gauge 
boson and a fermion is induced by the following terms in the Lagrangian, 

-£ = 4j dm Y [gi z P L + g$ z P R ] NZ, + h.c. , (3.3) 

- C = ^ DM 7 M [giwPL + 9*wPr] * + h.c. , (3.4) 

where Pl,r = |(1 =F 7 5 ) are the left and right-handed chirality projectors and the gfj R 
are coupling constants. These interaction terms yield the following decays widths for the 
two-body decays V'dm — > Z°N and ipDM ~^ W ± l T , 



and 
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respectively. From these decay widths we see that the couplings must be very strongly 
suppressed in order to yield dark matter lifetimes that exceed the age of the Universe. 

The decay of dark matter particles produces monoenergetic weak gauge bosons with 
an energy given by 



m DM 



and 



m DM A , M, 



Ew = ^ 1 + ^ , (3.8) 
respectively, where the charged lepton mass has been neglected. 

We show in Fig. [TJ in red lines, the constraints on the partial dark matter lifetime as 
a function of the dark matter mass in the case that the dark matter particle decays with 
100% branching ratio via ^dm — > W ± £ T . The long-dashed, short-dashed and dotted lines 
correspond, respectively, to the MIN, MED and MAX propagation parameters shown in 
Table [TJ The results for the decay ^dm _ > Z v are similar to the case -0dm — > W ± £ T , 
shown in Fig. [JJas orange lines, except in the region of dark matter masses which are close 
to the gauge boson mass. It follows from our calculation that the PAMELA measurements 
of the antiproton-to-proton fraction require, for dark matter particles which decay into 
weak gauge bosons, a corresponding partial lifetime longer than O(10 26 — 10 27 ) s at 
n^DM = 1 TeV, depending on the propagation model parameters when using the most 
conservative approach, and O(10 27 — 10 28 ) when using the more stringent and realistic 
approach which takes into account the background. 

Similarly, the two-body decay of a spin- 1/2 dark matter fermion with mass wdm 
into the Standard Model Higgs boson and a neutral particle is induced by the following 
interaction term in the Lagrangian, 

- £ = ^ M [\ L i!h P L + \% h P R ]Nh + h.c. , (3.9) 

where the \f/ R are coupling constants. This interaction yields a decay width given by 




+ ^Re(A^Aj*)l, (3.10) 



while the Higgs boson produced in the decay has an energy of 

^^(H^i (3.11) 

z V DM / 

In Fig. [TJ green line, we show the lower bound on the dark matter lifetime in the 
decay channel ^dm - > h°u as a function of the dark matter mass, where we assume that 
the Higgs boson is Standard Model-like and has a mass of m h o = 125 GeV. 
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Figure 1: Lower limits on the partial lifetime of a fermionic dark matter particle which 
decays into a lepton and a gauge or Higgs boson from the requirement of not exceeding 
the antiproton-to-proton fraction as measured by PAMELA at 95% C.L. The red line 
corresponds to the decay — > W ± £ T , the orange line to ^dm — > Z°iy and the green line 
to ^dm — > h°v, while the long-dashed, short-dashed and dotted lines correspond to the 
MIN, MED and MAX propagation model parameters, respectively (see text for details). 
Left panel: Constraints for primary antiprotons from dark matter decay + secondary 
antiprotons from cosmic-ray spallation. Right panel: Constraints for primary antiprotons 
from dark matter decay only. 

3.1.2 Three-Body Decays 

The energy spectrum of the antiprotons produced in three-body decays of dark matter 
particles contains some model dependence. In this paper we consider two possibilities, 
depending on whether the decay is mediated by a virtual charged scalar or a virtual 
charged vector particle, where we assume these particles to be much heavier than the 
dark matter mass. A generic interaction Lagrangian of the dark matter particle with a 
heavy charged scalar E which induces the decay ^dm — > qqN is 

£ cff = ~^DM[>^PL + X^P R ]q^-N[X^ N P L + Xf N P R ]qT^ + h.c. (3.12) 

In the limit m q <C ?tidm <C mj the partial decay width for the decay V'dm — > qqN is 
given by 

r(^ D M -+ qqN) = 128(2 ^ 3m 4 {C^FMIKm) + Cf F 2 (m 2 N / m 2 M )} . (3.13) 
The constants Cf, Cf are determined by the couplings as 

Cf = (|A^| 2 + |A* | 2 ) (|A^| 2 + \Xf N f) - V Re (A^A^A* A*) , (3.14) 
Cf ee 2 v Re [(A^A^) 2 + (A&A$) 2 1 . (3.15) 
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^dm(p) 




Figure 2: Tree-level diagrams contributing to the three-body decay ipr>u — ► qqN of 
fermionic Majorana dark matter, mediated by a heavy charged scalar S. Instead of the 
intermediate scalar £, the decay can also be mediated by a vector V. 



Here, r\ = rj^ DM rj N = ±1 depending on the CP eigenvalues of ^dm an d N. The kinematical 
functions, on the other hand, are given by 

Fi(x) = (1 - x 2 )(l + x 2 - 8x) - 12x 2 ln(x), (3.16) 
F 2 {x) = y/x[(l - x)(l + 10z + x 2 ) + 6x(l + x) ln(x)]. (3.17) 

In the hierarchical limit ttitvA^dm — > 0, the kinematical functions satisfy 

Fi(x)~l, F 2 (x)~v^, (3.18) 

whereas in the degenerate limit m N /m DM — » 1, one gets 

Fx(x)~ jj(l-x) 5 , F 2 (x)~^(l-x) 5 . (3.19) 

In the case of a vector interaction, on the other hand, we consider an effective La- 
grangian of the form 

= -^dm7 m [A^Pl + A« 9 V} - N-f [\ L qN P L + \* n Pr] q V* + h.c. (3.20) 

The decay width for the three-body decay mediated by a heavy charged vector particle 
reads 

r(^ D M -+qqN) = — g% 4 {C^K/r^M) + C 2 ^KK M )} , (3.21) 
32(27r) ,;l my 1 J 

where the functions Fi and F 2 were defined in Eqs. (I3.16p . (13 . 1 7[) and 

C\ = (|Aj/ + I A?/) (|A^| 2 + lAf^l 2 ) + 2 V Re (A^A&A&A*) , (3.22) 
C 2 y ee 2r/Re [(A^A^) 2 + = C 2 S . (3.23) 

10 



The energy spectrum of quarks and antiquarks produced in the three-body dark 
matter decay is fairly model dependent, since it depends both on the couplings which 
are involved in the decay as well as whether the mediating particle is a heavy charged 
scalar or a vector. More specifically, when the decay is mediated by a heavy charged 
scalar, the differential decay rate is: 

°^ _ 3 m DM -E 2 (-E max - E) 2 J L 2 R 2 L |2 I | \R |2\ v 

dE ~ 2(2vr)3 (m DM -2£)mf 1 " I'VI KI a 9 jvI " IVN x 



16£ 2 + 2E(E max - 9m DM ) - 3m DM (£ max - 2m DM ) 



3(m DM - 2£) 5 
+ r]Re (A^,A^r) + (A^,A^) 

- 2?7Re (A^A^A^A^) 



(3.24) 



m D M - 2£ 



where 



rn 



DM 



In the hierarchical limit, m^/miyyi — > 0, the normalized energy spectrum has a 
xin 
yields 



maximum and a minimum which is reached when A^ = A^, = A^y = and which 



1 dT\ n 2 .. / n f A 5 A 6 A 11 



with y = E/E ma _ x . On the other hand, in the degenerate limit mjv/mDM — > 1, the 
normalized energy spectrum is independent of the combination of couplings and reads 

If) ^30y 2 (l-y) 2 . (3.27) 



TdyJ 

The allowed ranges for the energy spectrum in the case of decays mediated by a heavy 
charged scalar in the hierarchical and degenerate limit are shown in Fig. [3j 

In the case that the decay is mediated by a heavy charged vector, the differential 
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decay rate is: 

dT 
dE 



6 m BM E 2 {E raax -Ef 



(|AiA 



L |2 
qip /s qN I 



(2tt) 3 (m DM - 2E)m A v 
l&E 2 + 2£(£ max - 9m DM ) - 3m DM (3 



, I \R \R \2\ 
\ A q%l) A qN\ ) 



2m 



DM J 



3(m DM - 2£) 2 



+ ?7Re (A^A^r) + (A^A^) 



,R \-R*\ 2 



+ 2 



(|AiA 



R |2 



,R \L 1 2> 



qip'\N\ \ A qip A qN 



AE 2 + 2E(E n 



m DM - 2£ 
+ 2?7Re (A^A^A^A^) 

3m D M) - 3m D M(-E max ^DM, 



3(m DM - 2Ef 



(3.28) 



In this case, the maximum and minimum of the energy spectrum is defined by the 
two curves 



i dr\ 



max, mm 



9y 2 x (max, min) \ ( 1 - -y ) , - ( 1 - - y 



(3.29) 



which correspond to the cases where the dark matter particle and the daughter neutral 
particle couple to quarks with the same chirality or with opposite chirality. On the other 
hand, in the degenerate limit we find: 



1 dT\ o, 



yf 



(3.30) 



which is identical to the result for the decay mediated by a heavy scalar particle. 

For the hierarchical limit, it is apparent from the figure that both spectra are very 
similar for E < 3/4E max , whereas at the endpoint of the spectrum the difference amounts 
at most to a factor of 2. The largest difference corresponds to the case where the dark 
matter particle and the heavy neutral daughter particle iV couple to quarks with different 

(or, alternatively, when 



chiralities, namely in the scenario where A^, 



and A^v 



A^, = and A^ = 0). In contrast, it can be checked from Eqs. (|3T24|l . (13T28I) that the 
decay spectra are identical when the dark matter particle and the heavy neutral particle 
N couple to quarks with the same chiralities. Since the spectra are fairly similar we only 
analyze the case with the decay mediated by a scalar and the chiral couplings, motivated 
by the case when the neutral daughter particle is a neutrino. 

We show in Fig. H] the constraints on the dark matter lifetime as a function of the 
mass in the case that the dark matter particle decays into a quark-antiquark pair and a 
neutrino. In the left-hand plot we require that the total flux, including the background 
antiproton flux calculated in [23], does not exceed the antiproton-to-proton ratio observed 
by PAMELA. In the right-hand plot we show conservative constraints from requiring 
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Figure 3: Energy spectra of the quarks produced in the three-body decays -0dm — > qqN 
mediated by a heavy scalar and a heavy vector, respectively, normalized to the total decay 
rate (blue: scalar, red: vector), when assuming that <C ttidm- The black dashed line 
shows the spectrum in the nearly degenerate case < tti-qmi which is identical for 
scalar and vector-mediated decays. 

that the contribution from dark matter decay alone does not exceed the antiproton-to- 
proton ratio. The red line corresponds to the constraints for the decay into light quarks, 
concretely the decay into a down-quark pair and a neutrino, ipDM — > ddu; the results for 
the decay into u and s quarks are very similar and are not displayed. On the other hand, 
the constraints from the decay into the heavy quarks c, b and t are shown, respectively, 
in green, orange and blue. Notice that the constraints from decays into heavy flavors are 
more stringent, because of the higher multiplicity of antiprotons in the fragmentation 
of heavy quarks compared to the decays into the lighter quark flavors. This effect is 
particularly important for decays into top-antitop pairs, however, this decay is only 
relevant for relatively large dark matter masses, > 2m t ~ 350 GeV, for which this 
decay is kinematically accessible. 

In Fig. [5] we show a comparison between the constraints resulting from the decay 
V'dm — > ddv with a massless fermion v and the decay ^dm — > ddN with a massive fermion 
N which is near-degenerate in mass with the dark matter particle, = 0.9 m^M- 
Instead of mass vs. lifetime, we show the maximum quark momentum vs lifetime times 
mass in this plot. Because the flux of antiprotons is proportional to PdmF/ttidm; the 
combination mDMT -1 can be constrained independent of the dark matter mass, allowing 
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for direct comparison between the massive and massless cases. The constraints are quite 
similar, although not completely identical due to the different energy spectra in the 
massless and massive cases. 




100 1000 10000 100 1000 10000 



"dm [GeV] m DM [GeV] 



Figure 4: Same as Fig. [TJ but for the decay of a fermionic dark matter particle into a 
quark-ant iquark pair and a neutrino. The red line corresponds to the decay -^dm — > ddu, 
green to t/'dm — > ecu, orange to ^dm — > bbu and blue to ipr>M — > tiu. 




[GeV] 



Figure 5: Direct comparison of the constraints on the nearly mass-degenerate three-body 
decay ipBM —* ddN with m N = 0.9?tidm and the decay ^dm - > ddv with a mass- 
less fermion v. The thin lines correspond to the massless case, whereas the thick lines 
correspond to the massive case. On the horizontal axis we plot the maximum quark mo- 
mentum, while we plot the combination mDMT" 1 on the vertical axis as explained in the 
text. 
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3.2 Decay of Spin-0 Dark Matter Particles 

In the case of a decaying scalar dark matter particle, we consider the following decay 
channels: 

• Decay into a pair of weak gauge bosons: 0dm — > Z°Z°, 0dm — > W + W~ . 

• Decay into a pair of Higgs bosons: 0dm -> h°h°. 

• Decay into a quark-ant iquark pair: 0dm — > QiQj- 

To calculate the antiproton fluxes at Earth we again follow the procedure outlined in 
subsection 13. 1[ 

We show in Fig. E] the resulting constraints for the scalar dark matter decay 
modes 0dm —> W ± W T , 0dm — > Z°Z° and 0dm — > h°h°, as red, orange and green 
lines respectively, where we again assume a Standard Model Higgs boson with a mass 
m h o = 125 GeV. 




1000 10000 1000 10000 

todm [GeV] WDM [GeV] 



Figure 6: Lower limits on the partial lifetime of a scalar dark matter particle which 
decays into two bosons from the requirement of not exceeding the antiproton-to-proton 
fraction as measured by PAMELA at 95% C.L. The red line corresponds to the decay 
0dm W ± W Zf , the orange line to 0dm — > Z°Z° and the green line to D m _ > h°h°, while 
the long-dashed, short-dashed and dotted lines to the MIN, MED and MAX propagation 
models, respectively (see text for details). Left panel: Constraints for primary antiprotons 
from dark matter decay + secondary antiprotons from cosmic-ray spallation. Right panel: 
Constraints for primary antiprotons from dark matter decay only. 

The constraints for decays into quark-antiquark pairs, on the other hand, are shown 
in Fig. [3 We show in red the constraints on the lifetime for the decay mode 0dm - > dd, 
in green for D m - > cc, in orange for D m - > bb and in blue for D m - > tt) the limits 
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for decays into the light quarks u and s are very similar to the limits for decays into d 
quarks and are not displayed. 
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Figure 7: Same as Fig. El but for the decay of a scalar dark matter particle into a quark- 
antiquark pair. The red line corresponds to the decay 0dm — > dd, green to 0dm — > cc, 
orange to 0dm _ > bb and blue to 0dm _ > ti. 

4 Constraints from Radiative Decays into Gamma 
Rays 

A dark matter particle that decays into electrically charged particles at tree-level gener- 
ically also possesses a decay mode producing monoenergetic photons that is induced 
by radiative corrections. As argued in [37], although the radiative decay mode is loop- 
suppressed, gamma-ray line searches can give rise to relevant constraints on the dark 
matter lifetime for leptonically decaying dark matter. Here, we extend this investigation 
to the case of hadronically decaying dark matter. More specifically, we analyze the con- 
straints on the parameter space of decaying dark matter scenarios where the dark matter 
particle has spin 1/2 and decays into a quark-antiquark pair and a neutral daughter par- 
ticle at tree level. In that case, a two-body decay producing a monoenergetic photon is 
induced at the one-loop level. The resulting constraint on the lifetime will thus be inde- 
pendent of the details of the propagation of the antiprotons from their production point 
to the Earth, which is the most important source of uncertainty in the constraints de- 
rived in the previous section. The constraints, however, will depend on the details of the 
underlying particle physics model, namely on whether the decay is mediated by a heavy 
charged scalar or by a heavy charged vector particle, as we will show below. 

We first consider the three-body decay of spin-1/2 dark matter particles into a quark- 
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antiquark pair and a neutral daughter particle, mediated by the exchange of a heavy 
charged scalar particle. The part of the Lagrangian which induces such decays is given in 
Eq. (13.12p . From this Lagrangian, the decay rate of the loop-induced process j^dm — > iN 
can be calculated in a straightforward manner. In the limit m N ,m q <C m DM <C wis the 
result reads 

3 



DM 



7 N) 



rn 



DM 



8vr (64tt 5 



rn 



N 



m 



DM 



rriBM 



X 



Qq { X qN X qip V X qN X q\) 



(4.1; 



where q = {u, d, c, s, b, t} runs over all quark flavors which interact with the dark matter 
particle and the neutral daughter particle, and Q q is the charge of the quark in the loop in 
units of the electron charge (Q q = 2/3 for up-type quarks and Q q = — 1/3 for down-type 
quarks) . 

Likewise, when the decay is mediated by a heavy charged vector, corresponding to 
the Lagrangian in Eq. ( 13. 2 1 f) . the decay rate of the loop-induced process ^dm — > lN 
reads 



9e 2 
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DM 
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rn 



DM 



87r(327r 2 ) z rri 



rjmN 

WDM 
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(4.2) 



where, as before, we have assumed the mass hierarchy m N ,m q <C m DM <C m 2 . 



4.1 Neutrinos as Neutral Daughter Particles 

Let us analyze first the case where the neutral daughter particle is a left-handed neutrino, 
implying A^y = and ~ 0. It then follows from Eq. (14.1 p that in order to obtain 
a sizeable decay rate into monoenergetic gamma rays, it is necessary that the dark 
matter particles also couple to quarks with left-handed chirality. Otherwise, the decay 
rate is proportional to the quark mass squared, which greatly reduces the intensity of 
the gamma-ray line unless the quark circulating in the loop is a top quark. If ^ 0, 
we find 

BR^dm — > 7^) — cm . qN g ^ ■ for intermediate charged scalar, 

BR(^dm — > 7^) — a ° m ■ -f=t ^It qN ,j Q to f° r intermediate charged vector. (4.3) 

87r l^q\ X qN X q^\ 
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We list in Table |2] the branching ratios for a series of interesting scenarios where the 
dark matter particle decays into just a single flavor of down-type or up-type quarks 
("0dm — > ddv and ^dm — > uuv respectively), democratically into all three quarks of the 
same type, either down type or up-type (■j/'dm — > QdQd 1 ^ an d V'dm — > QuQu^, respectively), 
or democratically into all quark flavors (^dm — > qqv). It is apparent from the table that 
in this scenario the relative branching ratio ranges between a era / (8n) ~ 3 x 10 -4 and 
108a cm /(87r) ~ 0.03, where decays mediated by virtual vector particles generally have 
a larger ratio of two-body decays to three-body decays than those mediated by virtual 
scalars. The PAMELA measurements on the antiproton to proton fraction constrains the 
inverse decay width of the three-body decays to be smaller than ~ 10 27 s, as shown in 
the previous section. Therefore, the inverse decay rate of the radiatively induced process 
ip — > '-/N will be in the range 10 29 — 10 31 s. The present observational upper bound on 
the inverse decay width into monoenergetic gamma rays is about 5 x 10 28 . . . 10 29 s for 
dark matter masses in the range 2 — 600 GeV. Furthermore, in the future the projected 
Cherenkov Telescope Array (CTA) [38] will constrain the inverse decay width in the TeV 
mass range also at the level 5 x 10 28 . . . 10 29 s from observations of the diffuse electron 
background, assuming an observational time of T = 1000 h and an energy resolution 
of 10% [37]. Therefore, for some concrete particle physics scenarios, the present and 
projected searches of the radiatively induced gamma-ray lines could provide constraints 
on those scenarios which are competitive with the ones stemming from the tree-level 
antiproton production. 

The interplay between gamma-ray constraints and antiproton constraints is illus- 
trated in Fig. [8] for the scenarios where the dark matter particle decays into down quarks 
and a neutrino (V'dm — ► ddv, top panel) or democratically into all three charge 2/3 quark 
flavors (ipDM qJivVi bottom panel), and which yield, respectively, the minimum and 
maximum ratio for the radiative two-body decays to the tree- level three-body decays (cf. 
Tab. [2]). In the plots, we show the constraints obtained in Section [3] together with the 
constraints on gamma-ray lines from two different analyses of the Fermi LAT data (dot- 
ted blue [39] , solid blue [ID]) as well as the prospected limits at the Cherenkov Telescope 
Array (cyan) [37]. In the plots, we use the upper limits on the inverse decay rate into 
monochromatic photons, and we scale them by the appropriate ratios for the radiative 
two-body decays to the tree-level three-body decays listed in Tab. [2j 

The limits published recently by the Fermi LAT collaboration [H] are very similar to 
the ones shown in Fig. [8] within the energy range of interest to us. Whether the feature in 
the gamma-ray spectrum near the Galactic center at E 1 ~ 130 GeV reported in p£2l H3] 
can be interpreted in terms of decaying dark matter remains to be seen [HI H5] , and we 
will not pursue this issue in further detail here. We just remark that the limits arising 
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27a om /(87r) 
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aorn / (87r) 


9a em /(87r) 


ip — y uuv 


4a em /(87r) 


36a cm / (87r) 




3a cm / (8vr) 


27a em /(87r) 




12a cm /(87r) 


108« cm /(87r) 


■0 — >■ ggi/ 


3a cm /(167r) 


27a cm /(167r) 



Table 2: Ratio r(^> — > 7^)/ J^/IX^ — >■ //^) of the one-loop decay rate into monoener- 
getic photons to the tree-level decay rate into fermions, for the case of scalar and vector 
mediated decays (neglecting fermion masses, and assuming purely chiral couplings). The 
decay channel ip — y £ + £~u corresponds to a dark matter particle which decays into 
e + e~u, fi + fi~u and t + t~i> with equal rates (flavor-democratic decay). Similarly, the last 
three lines correspond to flavor-democratic decay into the down-type quarks, the up-type 
quarks, or all quarks, respectively (assuming that ttidm 2m 4 in the latter two cases). 
The ratios are in the permille or percent range (a em = e 2 /(47r) ps 1/137). 

from hadronic decay channels discussed here constitute an important consistency check 
for any model aiming to explain the gamma-ray feature. 

4.2 Degenerate Scenario 

In concrete particle physics models, the neutral daughter particle produced in the three- 
body decay might not be massless, but could instead have a sizeable mass, possibly 
comparable to the dark matter mass. If this is the case, the branching ratio of the 
radiative decay into photons could be significantly enhanced by kinematic effects, as 
discussed in [37]. 

In the case of chiral couplings, for instance when the dark matter particle and the 
neutral daughter particle both couple to quarks with left-handed chirality, the branching 
ratio into monoenergetic gamma rays reads 

oo// 45a em / m 2 N \ ~ 2 f rim N \ 2 
BR ^dm -> 7^ ^ 1 r~ 1 

8?r V m E)M/ V m DM / 





Y^q Qq { A qN A qi) 


- 2 


\^\ A qip\ 


2 |Ay 2 + 77Re 


(\L \L*\ 2 
\ A qip A qN) 


} 



(4.4) 
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Figure 8: Comparison of constraints from cosmic-ray antiprotons and radiatively gen- 
erated gamma-ray lines for the tree-level decay into a pair of down-type quarks and a 
neutrino (^dm — >* ddis, top panel), as well as for the democratic decay into charge 2/3 
quark flavors (^dm — > QuQu^, bottom panel). The thick black lines are the limits on 
the decay width from antiprotons including the background from cosmic-ray spallation, 
while the thin lines correspond to the limits resulting from dark matter decay only. On 
the other hand, the blue lines represent the limits from gamma-ray lines from the Fermi 
LAT data [39l HQ], while the cyan lines represent the projected limits for the Cherenkov 
Telescope Array [37] . Left panels: Intermediate scalar. Right panels: Intermediate vector. 
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when the decay is mediated by a heavy charged scalar, and 



BR(^ DM -)• -yN) 




(4.5) 



E 9 { 2 l\".l 2 l A ^l 2 + ^ Re [W^) 2 ]} 



when the decay is mediated by a heavy charged vectorjj In both cases, one finds a 
significant enhancement of the decay width into monoenergetic gamma rays as mjy -> 
777. dm when the dark matter particle and the daughter particle have opposite CP parities, 
77 = — 1. In this case the branching ratio has the proportionality 



which can lead to a significant enhancement of the radiative decay as — > ^dm- 

In Fig. [9] we show the constraints from gamma-ray and charged cosmic-ray obser- 
vations for the case that the dark matter particle ipDM and the neutral fermion N are 
nearly mass-degenerate, = 0.9 777dm, and have opposite CP parities, 77 = —1. Due to 
the enhancement of the radiative decay mode in this case, the lower limits coming from 
the loop-induced decay into monochromatic photons are comparable to the lower limits 
obtained from the primary antiproton flux produced in the tree-level decays for scalar- 
mediated decays and for MIN or MED propagation parameters. For vector-mediated 
decays, the gamma line limits are even stronger than the antiproton limits obtained for 
the MAX propagation parameters. 

5 Constraints from Radiative Decays into Weak 
Gauge Bosons 

It has been proposed that the observation of cosmic ray anomalies in positron and elec- 
tron fluxes reported by the PAMELA and Fermi satellite experiments could be due to 
dark matter that decays or annihilates primarily into leptonic channels jlHl H7] (see, 
e.g., Ref. [6] for a recent analysis of this scenario). Nevertheless, even leptophilic dark 
matter will lead to a certain amount of antiproton production induced by electroweak 

^^Note that we have used the asymptotic approximation Eq. (|3.19[) for the kinematical functions Fi(x) 
for obtaining these expressions, i.e. they are valid for tojv — wdm- We have checked that, using the 
full Fi(x), the resulting branching ratio smoothly interpolates between the expressions given here for 
raw — mdm and those for mjy — > derived in the previous subsection. 




(4.6) 



21 




1000 10000 1000 10000 



m D M [GeV] bidm [GeV] 



Figure 9: Same as Fig. El top panels, but for the near-degenerate case = 0.9 moM- 
Comparison of constraints from cosmic-ray antiprotons and radiatively generated 
gamma-ray lines for the tree-level decay into down-type quarks ("0dm — * ddN) in the 
degenerate case where = 0.9 mDM- The corresponding gamma-ray constraints are 
stretched in the moM-direction because the available gamma-ray energy is lower if N is 
massive. Furthermore, the constraints are scaled in the T _1 direction due to the kine- 
matic enhancement resulting from the opposite CP parities of V'dm and N, as discussed 
in the text. 

bremsstrahlung, as has been argued, e.g., in Refs. [381 139"]. We consider an additional 
source of antiproton production resulting from radiative decays into weak gauge bosons. 
We first consider the case of a fermionic dark matter particle and then the case of scalar 
dark matter. 

5.1 The Decay ^ DM -> Z°N 

We consider a fermionic dark matter particle that decays into a pair of leptons and a 
neutral particle, -0dm — > £ + £~N, at tree level. At one-loop, the decay -0dm — > lN is 
induced, allowing for constraints from gamma-ray searches [37]. In addition, the decay 
■0dm Z°N into a Z boson is generically induced, which gives rise to antiproton pro- 
duction. Depending on the model, the decay mode ipr,M — > W ± £ T can also be induced. 
In order to study the radiative decays into weak gauge bosons, it is necessary to specify 
the S , 6 r (2)i,-structure of the model. Here we investigate a scenario where the dark matter 
particle V'dm and its neutral decay product iV are both Standard Model singlets, and 
consider the effective interaction Lagrangian 

C cS = -AL^dmLE 1 - XfMulR^ - \ L m NLtf - \f N N£ R ^ + h.c. , (5.1) 
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where L = \ vz iij is the left-handed lepton doublet and £r is the right-handed charged 
lepton field. The heavy scalar fields £ and a mediate the decay, and have the same 
quantum numbers as L and £r, respectively. In a supersymmetric context, they can 
be identified with the left- and right-handed slepton fields. The effective Lagrangian 
leads to an interaction of the form described by Eq. (13.121) . except that two mediating 
fields, namely the charged component of £ as well as a, contribute. In addition a similar 
interaction involving neutrinos and the neutral component of E is present. For simplicity, 
we neglect mixing of the mediating fields, i.e. we assume that £ and a coincide with mass 
eigenstates, and we also assume that the vacuum expectation value = 0. 



Figure 10: Diagrams contributing at one loop to the radiative two-body decay ^dm — > 
Z°N, induced by a charged scalar E (top row) and a vector particle V (bottom row), 
respectively. There are two additional diagrams in each case which differ only by the 
direction of the charge flow. 

The decay rates for the tree-level decays -0dm — > £ + £~N and ^dm — > vvN as well 
as the loop- induced decay t/>dm - > iN can be adapted from Ref. [31]. In addition, in 
this scenario the decay -0dm - > Z°N occurs at the one loop level, leading to antiproton 
production from Z boson decay. The corresponding Feynman diagrams are shown in the 



^dm(p) 
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upper row of Fig. [TU1 We find for the decay rate in the limit mAr, <C ttidm, Mz 
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with ms being of the order of the intermediate particle masses^] 

For example, for purely left-chiral couplings the branching ratios for loop-suppressed 
decays are given in this scenario by 
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(5.4) 

where the three cases refer to the possible mass orderings of the charged and neutral 
components of the mediating scalar fields. For comparison, the branching ratio of the 
electromagnetic decay mode for these cases is given by 



BR(^ DM jN) 



3a em I 52e A w A 



tip\ 



87T EMn A ^\ 2 
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4 

— r 2 - < 1 m s± > m So 



(5.5) 



As expected, both decay channels differ mainly due to the different coupling strength, 
as well as the different kinematics arising from the finite mass of the Z boson. When 
m^ ± = ms , the dark matter particle decays with equal rate into leptons and neutrinos 
at tree level. Therefore, the branching ratio for the decay into photons is half as large 
compared to the case ms ± <C ms . In the case m^ ± ^> ms the electromagnetic channel 
is relatively suppressed, while the radiative decay producing Z bosons is also possible 
via the vertex diagrams involving electrically neutral particles in the loop. 

The diagrams for the case when the particle mediating the decay is a vector boson 
are shown in the lower row of Fig. [TQ1 The resulting branching ratios can be obtained by 



2 In general the masses appearing in the logarithms in F% can be different for each mediator. We 
assume here for simplicity that the masses are of similar size. Alternatively, if the rate is dominated by 
one intermediate particle, then ms corresponds to the mass of this field 
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Table 3: Benchmark scenarios. In the first three cases, the three-body decay produces 
electrons only. In the last four cases, the three-body decay produces muons only. 



replacing in Eq. (I5.2I) m s — > m v , 3a em /(87r) — > 27a cm /(87r) and Ff — > F%, where 

Fj(i) = 1 + ! (l - | (Hm 2 v /Ml) + 3) + + 2x) \3 + Hm 2 v /M 2 z ) - m| 2 ^) . 

(5.6) 

As for the case of the loop-induced decay ip — > jN, the rate for the process ip — > ZN 
can be enhanced in the degenerate limit. This is illustrated in Fig. [TT1 where we plot the 
constraints on mDMr -1 against the energy of the Z boson. In the figure we also display the 
positions of the benchmark points from Table [3] for comparison. We see that depending on 
the ratio between 7Bdm and m N , the radiative decays can vary in their relative intensity 
to lie either in the allowed or in the excluded region of the parameter space. Benchmark 
points 1 and 4, which correspond to the highest degeneracy and both have r] — — 1, 
can be ruled out based on the primary antiproton flux produced from the loop-induced 
decay into Z bosons, when assuming MAX or MED propagation parameters, respectively. 
Points 3 and 6 also have an enhanced radiative decay, but are still in the allowed region. 
For comparison, points 2 and 5 (for which 77 = +1, so that there is no enhancement 
of the one-loop decay) are also shown. Note that the limits arising from the one-loop 
decay into monoenergetic gamma rays have been analyzed for the same set of benchmark 
points in [37]. There, it was found that the photon flux predicted for benchmark points 
1 and 4 is in conflict with Fermi LAT limits from the Milky Way halo, and MAGIC 
observations [50] of the Perseus cluster, respectively. The photon flux of benchmark point 
6 can be tested in the future by the CTA. Thus, we conclude that for dark matter that 
decays "leptophilically" at tree level, the limits from monochromatic gamma rays and 
from antiprotons produced in the one-loop decay channels into photons and Z bosons, 
respectively, are competitive, although the former ones are slightly more restrictive in 
the case of MED and, especially, for MIN propagation parameters. 



25 





■ B 5 


■ .B 2 






■ B 6 ~ 








7 -B 4 " 





100 1000 

E z [GeV] 



Figure 11: Constraints on the decay "0dm Z°N with > 0. We show the benchmark 
points from Table [3] for comparison. 



5.2 The Decay DM -> 7^° 

Analogously to the fermionic case, we also consider a scalar dark matter particle that 
decays leptonically via its tree level coupling to charged leptons, 0dm — > £ + £~. At the 
one-loop level, the dark matter can also decay into a pair of gauge bosons. The rate for 
the loop-induced decay into a pair of photons, 0dm — > 77, has been discussed in the 
context of gamma-ray line searches in Ref. [37]. Antiprotons are produced in the loop- 
induced decays involving the massive gauge bosons. The decay rate for the loop-induced 
process 0dm _ > lZ° can be adapted, e.g., from the results in Refs. [511 [52]. For the 
branching ratio, we find in the limit rri£ <C ttidm, Mz 



pp/i ^„ /7 on 4aLm, 2 (l-x) (3s 2 w - c 2 w 

hJK(0DM -> 7^ ) — 5 2 2 — 2 



f2m DM S W C W 



3/1 x 
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IO-^TT^Tt) ^ (1-5^) , (5-7) 



xl06 MeV/ \ ™dm / V m bM 

where x = M^/m 2 :)M and xi = m^ DM /m 2 . This radiative decay mode is suppressed 
because of the necessary helicity flip of the lepton in the loop. Analogously, the decays 
0dm Z°Z° and 0dm — > W + W~ are in principle allowed. We do not examine these 
decay modes in detail since they are suppressed in the same manner. 
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6 Conclusions 

We have presented a general model-independent analysis of the cosmic-ray antiproton 
constraints on the partial dark matter lifetime in hadronic decay modes for the lowest- 
order decays allowed by Lorentz and gauge invariance. In particular, we have investigated 
the possible two-body decay modes for fermionic and scalar dark matter particles, as well 
as three-body decays into a pair of quarks and a neutral particle (for example, a neu- 
trino) for Majorana dark matter particles. For the three-body decay we have paid special 
attention to the resulting energy spectra, which can differ slightly depending on whether 
the decay is mediated by a scalar or a vector particle, and whether the neutral particle 
produced in the final state is nearly massless or nearly degenerate with the dark matter 
particle. Employing a semi-analytical model of cosmic-ray propagation we have scanned 
the mass-lifetime parameter space over several orders of magnitude, and have derived 
lower limits on the lifetime for three representative sets of propagation parameters from 
the PAMELA measurement of the antiproton-to-proton ratio. We have also examined 
kinematic enhancement of radiative decays as well as the production of monochromatic 
gamma rays. Radiative decays into photons induced at next-to-leading order by tree- 
level decay into the electrically charged quarks provide an interesting complementarity 
to the antiproton constraints. Finally, we have discussed antiproton constraints arising 
from loop-induced decays into weak gauge bosons for leptophilic dark matter. 
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Note Added 

During the completion of this work an independent preprint [53] appeared in which 
three-body dark matter decays and antiproton constraints are discussed. 
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